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Abstract. We have performed VLT spectroscopy of an almost complete sample of 18 early-type galaxies with 0.88 < z < 1.3 
plus two at z = 0.67, selected from the K20 survey, and derived the velocity dispersion for 15+2 of them. By combining 
these data with HST and VLT images, we study the Fundamental Plane (FP), the Faber-Jackson and the Kormendy relations 
at z~l, and compare them with the local ones. The FP at z~l keeps a remarkably small scatter, and shows both an offset and a 
rotation, which we interpret in terms of evolution of the mass-to-light ratio, and possibly of the size. We show evidence that 
the evolution rate depends on galaxy mass, being faster for less massive galaxies. We discuss the possible factors driving the 
evolution of spheroids and compare our results with the predictions of the hierachical models of galaxy formation. 
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1. Introduction 

The K20 survey (Cimatti et al. 2002c ) with its selection in the 
near infrared and its very high spectroscopic redshift complete- 
ness (92%) provides an unprecedented tool for studying the his- 
tory of galaxy mass assembly in the Universe, and in particu- 
lar the formation of massive galaxies. In the currently popular 
ACDM scenario, massive galaxies are the product of rather re- 
cent hierarchical merging of pre-existing disk galaxies taking 
place largely at z < 1.5 and with moderate star formation rates 
(e.g. Kauffmann et al. 119931 and Cole et al. I2000> . Therefore 
in such hierachical merging scenario fully assembled massive 
galaxies with Ai > 10 n vV(o and with evolved stellar popula- 
tions should be very rare at z > 1 . The alternative possibility 
is that massive systems formed at much higher redshifts (e.g. 
z > 3), through a short and intense period of star formation, 
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followed by passive evolution of the stellar population. Such a 
possibility is supported by the properties of local and interme- 
diate redshift spheroids (Renzini 1999 1, and by the existence 
of old and passive galaxies at z~1h-2 (e.g. Cimatti et al. I2002al 
Cimatti et al. 120041 McCarthy et al. 120041 Saracco et al. l2005l . 
The controversy can then be solved by finding a number of 
evolved massive galaxies at z > 1 . The problem so far has been 
the reliable determination of the mass, since the stellar masses 
estimated using the mass to infrared luminosity ratio or the fit to 
the multicolor spectral energy distribution are model dependent 
(IMF) and subject to various degeneracies (age - metallicity - 
extinction). 

The advent of 8-10m class telescopes has allowed dynam- 
ical masses to be estimated more reliably by kinematic stud- 
ies of galaxies at z~l. Various groups (e.g. van Dokkum & 
Stan fordl20Q3l Treu et al. l2005bl H olden et al. l20Q5l Gebhardt 
et al. 2003 van der Wei et al. 2005 and references therein) have 
conducted kinematic studies of early-type galaxies, selected by 
optical morphology, both in the cluster and in the field environ- 
ment up to z~ 1.3. They have used the Fundamental Plane (FP, 
Djorgovski & Davis 1987 and Dressier et al. 1987) to anal- 
yse the evolution of the Ai/L ratio, generally finding that the 
most massive spheroids must have formed at rather high red- 
shift (z~3). In order to understand the formation and evolution 
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of early-type galaxies over a broader mass range it is neces- 
sary to extend these studies with different selection criteria and 
better completeness. 

We address this problem by making spectroscopic obser- 
vations of a sample of high redshift galaxies selected from the 
K20 survey, and by exploiting the large collecting area of the 
ESO Very Large Telescope (VLT), enhanced by the exception- 
ally good red sensitivity of FORS2 on the VLT. The selection 
of the K20 survey, being in the K s band, which approximately 
corresponds to the rest frame /-band at the redshift of our sam- 
ple, is much more sensitive to the stellar mass of the galaxy and 
much less dependent on the possible presence of star forming 
activity than the usual selection in the observed optical bands. 
Therefore it is suited to study the evolution of massive galaxies 
in a less biased way. Furthermore the spectroscopic classifica- 
tion of the K20 survey assigns a galaxy to the early-type class 
based on the presence of old stars, definitely a much more sta- 
ble property than morphology and therefore less likely to ex- 
clude possible progenitors of today's spheroidal galaxies. In 
fact the only progenitors that we might miss are those where 
the clear spectroscopic signature of old stars would be over- 
whelmed by a major starburst, involving a considerable frac- 
tion of the galaxy mass. Our sample should therefore be less 
affected by the so called "progenitor bias" (e.g. van Dokkum & 
Franx 2001 1. Finally, the fact that the K20 survey and our study 
cover two independent fields, well separated in the sky, helps 
in reducing the effects of cosmic variance. 

We describe here the results of our kinematic study on the 
early-type galaxies of the sample selected from the K20 sur- 
vey, while a parallel paper (Vernet et al., in preparation) reports 
on those for disk galaxies. We assume a flat Universe with 
Q m = 0.3, Qa = 0.7, and Ho = 70kms~ 1 Mpc~ 1 , and we use 
magnitudes based on the Vega system. 

2. Sample selection and observations 

The galaxies to be observed have been selected from the 
K20 survey, which contains 545 objects selected in the K s 
band, with A' J <20.0, in two areas of the sky, the CDFS and 
0055-269 fields covering 32.2 and 19.8 arcmin 2 , respectively. 
We have selected galaxies with spectroscopic redshift between 
0.88 and 1.3 and with early-type spectra, i.e. classified as 
class 1.0 (early-type galaxy without emission lines) and class 
1.5 (early-type galaxy with emission lines) in the original 
K20 classification, based on a visual inspection of the spectra 
(Cimatti et al. 2002c I. The limits of the redshift range, partic- 
ularly the lower one, were set by the need of covering with a 
single spectroscopic setting the rest-frame wavelength interval 
between 3700A and 4500A, which contains the most useful 
absorption lines for the measurement of the velocity dispersion 
and the [Oil] 3727 emission line doublet. 

Of the 14 and 10 galaxies selected with our criteria in the 
CDFS and 0055-269 fields respectively, we have observed 
9 and 9 objects (see Table [TJ, i.e. 75%. The remaining ob- 
jects have been excluded for purely practical reasons, i.e. they 
were either outside of the mask prepared for multi-object spec- 
troscopy or their position was conflicting with that of other 
objects included in the mask. A better completeness has been 



achieved on the 0055-269 field because this, being smaller, is 
easier to cover with a single FORS2 mask. The result of our 
good completeness is that the observed galaxies span across 
the whole range of luminosities available for K20 early-type 
galaxies and give an unbiased view of the properties of early- 
type galaxies in the selected redshift and magnitude ranges. In 
order to check the properties of galaxies at lower redshift, we 
have included in the observed sample two galaxies at z = 0.67, 
one for each field. Recently a more quantitative classification 
of the K20 spectra has become available (Mignoli et al. 2005 1, 
with some differences with respect to the original visual clas- 
sification. The only practical effect on our sample is that the 
galaxy identified as CDFS_00547 (which is class 1.5 in Cimatti 
et al. 2002c I, with the recent reclassification would have been 
excluded from the sample, since it is not classified as an early- 
type galaxy, because its 4000A break is too small, although 
the spectral energy distribution is red (Mignoli et al. 2005 1. Of 
the 20 observed galaxies (see Tabled 1 1 qualify as extremely 
red objects (ERO), according to the R — K s > 5.0 definition 
(Elston, Rieke & Rieke 1988), a confirmation of the power of 
this colour selection in identifying evolved galaxies at interme- 
diate redshift. 

The selected redshift range excludes the two peaks in the 
redshift distribution around z~0.67 and z~0.74 in the K20 sur- 
vey (Cimatti et al. 2002b). Therefore our sample is most likely 
made of field galaxies, except for the two galaxies at z = 0.67, 
outside our main redshift range, which probably belong to the 
structures present at that redshift in both fields. The observed 
sample of galaxies is listed in Tabled The reader is referred to 
the recent release of the K20 survey 1 for the object identifica- 
tion and for additional information on the individual galaxies. 

The spectroscopic observations to obtain the velocity dis- 
persion have been performed between 31 October and 2 
November 2002 with FORS2 on UT4 of the ESO VLT at 
Paranal. We have used the grism 600z with slits 1.0 arcsec 
wide, which gives a resolution A/AA = 1400 and covers the 
wavelength range 7400-10700A. The rest-frame wavelength 
range for each individual galaxy, although varying according 
to its redshift and its position in the field, includes always the 
H and K [Call] lines and the G band, except for one of the 
two lower redshift galaxies (q0055_00169), for which the spec- 
trum starts around 4000A. We used two multi-object masks, 
one for each field. The exposure time was 7.5 and 6.0 hours, 
and the average seeing was 0.7 and 0.625 arcsec for the CDFS 
and 0055-269 field respectively. For calibration we have ob- 
served the spectrophotometric standard stars Hiltner 600 and 
HD49798. 

For the morphological and photometric analysis we have 
used the HST+ACS 5 epochs mosaicked stacks F850LP 
images from the Great Observatories Origin Deep Survey 
(GOODS, Giavalisco et al. 2004) and images from the ESO 
Imaging Survey (EIS, NTT+SUSI2 U,B, V,Rmd NTT+SOFI 
J, K s ) for the CDFS field, and images obtained with the 
NTT+SUSI2 (U,B,VJ), the NTT+SOFI (J,K S ) and the 
VLT+FORS1 (R,z) for the 0055-269 field. 



1 See http://www.arcetri.astro.it/~k20/spe_release_dec04/index.html 
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Table 1. The sample of observed early-type galaxies from the K20 survey 



K20 name 


z 


R 


R - K s 


n 


K v 


Re 

kpc 


<H B ) e 
mag/arcsec 2 


O" 

km/s 


EW[OII] 
A 


M B 


JTT, 

u 


CDFS.00060 


1.188 


24.89 


5.56 


















CDFS.00369 


0.8930 


23.30 


4.14 


1.85±0.39 


7.37+0.18 


0.84+0.05 


17.78+0.15 


119+21 


9.0+1.1 


-20.41 


10.31+0.15 


CDFS.00467 


0.8956 


22.90 


4.63 


3.45±0.39 


5.72+0.25 


1.52+0.10 


18.40+0.15 


140+18 


4.3+0.8 


-21.09 


10.58+0.11 


CDFS.00468 


1.019 


24.86 


5.03 


















CDFS.00532 


1.2115 


23.92 


5.22 


2.04±0.33 


6.82+0.20 


1.41+0.10 


17.43+0.16 


260+30 


2.5+0.5 


-21.89 


11.16+0.10 


CDFS.00547 


1.2243 


23.68 


5.07 


1.90+0.31 


7.50+0.14 


0.74+0.04 


16.15+0.14 


256+28 


2.3+0.4 


-21.76 


10.93+0.10 


CDFS.00571 


0.9551 


22.33 


4.79 


5.00±0.15 


3.83+0.02 


6.81+0.48 


20.65+0.16 


182+21 


3.7+0.5 


-22.08 


11.23+0.10 


CDFS.00590 


1.2235 


24.35 


5.28 


4.98±0.25 


4.03+0.12 


4.09+0.52 


20.54+0.29 


119+49 


4.1+1.1 


-21.09 


10.68+0.36 


CDFS.00633 


1.0963 


22.43 


5.49 


4.78±0.21 


4.00+0.12 


6.65+0.12 


19.97+0.05 


260+23 


< 1.4 


-22.71 


11.55+0.08 


CDFS.00354 


0.6672 


22.01 


4.07 


3.58+0.28 


5.20+0.20 


3.08+0.12 


20.39+0.09 


99+19 


2.6+0.6* 


-20.61 


10.52+0.17 


q0055_00028 


1.0524 


24.19 


5.09 


0.8±0.8 


7.17+0.48 


1.17+0.65 


18.08+1.22 


227+30 


< 3.6 


-20.82 


10.99+0.26 


q0055_00068 


1.1043 


23.53 


5.18 


4.3±1.0 


4.55+0.73 


3.60+1.00 


19.75+0.61 


84+31 


2.3+0.7 


-21.60 


10.37+0.34 


q0055_00114 


0.8891 


23.48 


4.33 


2.0+2.0 


6.71+1.78 


2.43+0.74 


19.99+0.66 


52+22 


< 2.1 


-20.50 


9.98+0.38 


q0055_00123 


0.9270 


22.96 


4.84 


3.8+1.0 


5.69+0.67 


1.23+0.40 


17.62+0.71 


186+19 


1.6+0.5 


-21.39 


10.74+0.16 


q0055_00247 


0.931 


24.43 


4.54 












<2.8 






q0055.00295 


1.1680 


23.75 


5.59 


4.8+0.3 


4.35+0.20 


2.67+0.70 


18.93+0.57 


170+28 


< 1.4 


-21.77 


10.85+0.18 


q0055.00318 


0.8962 


23.87 


5.07 


4.5+1.1 


5.17+0.70 


1.26+0.57 


18.87+0.99 


105+30 


< 0.8* 


-20.20 


10.21+0.31 


q0055_00331 


0.8980 


22.95 


5.09 


4.9+0.3 


4.15+0.19 


3.83+0.77 


20.08+0.44 


216+46 


< 0.8* 


-21.40 


11.18+0.19 


q0055_00338 


0.9324 


23.12 


4.96 


2.8+0.6 


5.93+0.50 


2.41+0.43 


19.23+0.39 


113+19 


<2.5 


-21.24 


10.59+0.16 


q0055_00169 


0.6699 


22.70 


4.23 


1.3+0.8 


7.15+0.62 


1.64+0.36 


19.50+0.47 


90+19 


< 0.9* 


-20.14 


10.32+0.20 



* for [OIII]5007. 



3. Data Analysis 

3.1. Morphology 

In order to obtain the effective radius R e , we have fitted a Sersic 
profile (logl(R) cc -R l l", Sersic fT%8l using GIM2D, a fitting 
algorithm for parametric two-dimensional models of surface 
brightness distribution (Simard et al. 1998). GIM2D performs 
a profile fit by deconvolving the data with the point spread 
function. We model PSFs with analytic functions from visually 
selected stars in the surrounding (30" x 30") region of each 
galaxy. We model a different PSF for each region in order to 
take account of PSF variations with the position in the field. 
A 2D radial gaussian function has been fitted simultaneously 
on tens of selected stars around the galaxies of our sample and 
outputs have been stacked together to result into a single PSF 
image for each region. A more detailed description of our ap- 
proach in modellig GOODS/ACS galaxy morphologies in the 
1.0 < z < 1.5 range will be given in Rettura et al. (2005, in 
preparation). 

The results of the bidimensional fit is the semimajor axis a e 
of the projected elliptical isophote containing half of the total 
light, the axis ratio b/a and the Sersic index n, which we have 
left as free parameters. The effective radius is computed from 
R e - a e yfbja and is given in Tabled 

Magnitudes have been measured on the ground-based im- 
ages using the SExtractor (Bertin & Arnouts 1996) BEST es- 
timator. The average surface brightness within the effective ra- 
dius (in mag/arcsec 2 ) is obtained from the absolute magnitude 
M: {y) e - M + 5logR e + 38.567, with R e in kiloparsec. 

In order to obtain the morphological parameters in the rest 
frame B band we have used the HST+ACS images taken with 



the F850LP filter and the VLT+FORS1 images taken with the 
Gunn z filter, since these are very close to the B band at the 
redshift of our galaxies. In any case a small K-correction has 
been applied following the prescriptions of Hogg et al. ( 2002 1 
and Blanton et al. ( 2003 ) and using the excellent information 
which we have on the spectral energy distribution. The surface 
brightness has been corrected for the cosmological (1 +z) 4 dim- 
ming. We have also corrected for the small Galactic extinction 
as obtained from Schlegel et al. ( 1998 1. 

For the 0055-269 field we lack HST images, so we have to 
rely on very good quality VLT+FORS1 images (seeing 0.62 
arcsec). Obviously the quality of the parameters obtained is 
worse than for those on the CDFS field. In order to check 
that the results are still acceptable and that the ground-based 
parameters are not biased, we have repeated the morpholog- 
ical analysis of the CDFS galaxies on a VLT+FORS1 image 
equivalent to that used for the 0055-269 field and found that 
the parameters obtained are consistent with those based on the 
HST+ACS images, although with larger errors (see Fig.Q. 

In order to check our consistency with the results obtained 
by others fitting the de Vaucouleurs law, we have repeated our 
fits by fixing the Sersic index to 4, obtaining sizes which are 
within 40% of those obtained with a variable index, some big- 
ger, some smaller, but without a noticeable dependence on the 
n value. The difference of the mean sizes derived with the two 
different fits is less than 5%. 

3.2. Spectroscopy 

The spectra have been reduced using the IRAF software pack- 
age in the usual way, including atmospheric extinction correc- 
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Fig. 1. Comparison of effective radii obtained for the CDFS 
field on the HST+ACS image and on the VLT+FORS1 one. 

tion and flux calibration. For each galaxy we have extracted a 
one-dimensional spectrum by coadding a section centred on 
the nucleus and one arcsecond long along the slit. This results 
in a spectral aperture of lxl arcsec 2 . The extracted spectra are 
shown in Fig.|5] 

The most noisy - i.e. sky contaminated - sections of the 
spectra have been masked out and the velocity dispersions have 
been obtained with the xcor task in IRAF, which uses the cross 
correlation technique of Tonry and Davis ( 1979 1. The correla- 
tion has been performed with template spectra of stars from the 
STELIB library (Le Borgne et al. I2003i . with spectral types 
ranging from F6 to K2, by matching their spectral resolution to 
that measured on our spectra. We have made tests with spectra 
of template stars artificially smoothed to a known velocity dis- 
persion, and with realistic noise added. These tests confirm that 
the results are not biased and that the quoted errors are realis- 
tic. The measured velocity dispersion has been normalized to 
a circular aperture with a diameter of 1.19/z -1 kpc, equivalent 
to 3.4 arcsec at the distance of the Coma cluster, following the 
prescriptions of J0rgensen et al. ( 1995bJ, as commonly used in 
FP studies. The changes on the velocity dispersion produced 
by these normalization factors are in any case small, i.e. 6-7%. 
The normalized velocity dispersion is listed in Tabled except 
for CDFS.00060, CDFS.00468 and q0055D0247, for which 
we could not obtain a reliable value because of the low S/N 
ratio. These three galaxies are the faintest in the R band. The 
average redshift of the galaxies with measured velocity disper- 
sion is z = 1.071 for the CDFS field and z = 0.983 for the 
0055-269 one, excluding the two galaxies at z = 0.67. 

Since we have noticed that several galaxies show a clear, 
although weak, [Oil] 3727 line emission, we have measured its 
rest-frame equivalent width (see Table 0. For the four galax- 
ies, for which the [Oil] line is outside the observed wavelength 



Fig. 4. The correlation of the Sersic index with the luminosity. 
Triangles are the Coma cluster data of Gutierrez et al. (|2004 1, 
black circles and red squares are our z~ 1 galaxies in the CDFS 
and 0055 - 269 fields respectively and empty symbols are for 
galaxies with z<0.88, while filled ones are for our redshift 
range 0.88<z<1.3. 

range, we have measured the [OIII]5007 line emission instead, 
as noted in the Table. 

4. Results 

4.1. Morphology 

If one expects that all early-type galaxies should have a de 
Vaucouleurs profile, then the Sersic indices that we obtain (see 
Tabled appear rather small, and there are three galaxies in the 
CDFS and three in the 0055-269 field with n < 2.5. However 
also in the local Universe several early-type galaxies are found 
to have a small n (see Fig.0}. All the galaxies in our sample do 
have a prominent spheroidal component, although a disk is vis- 
ible in a few of them (see Figure[3}, and all those in the CDFS 
field have been classified as E/S0 in the morphological analy- 
sis by Cassata et al. J2005I I. except for CDFS_00571, which is 
classified as Sa. 

Moreover we find that the Sersic index correlates with the 
B band absolute magnitude (see Figure^}, as found at low red- 
shift (e.g. Gutierrez et al. 2004 1, but with a shift to brighter 
luminosities for the same values of n. 

The distribution of the effective radii for our z~l galaxies 
is shifted towards lower sizes with respect to low redshift sam- 
ples (see Fig.[5]and Fig. II 0i. This could be due to a selection 
effect induced by the cosmological (1 + z) 4 dimming in sur- 
face brightness, in the sense that because of the dimming we 
select higher surface brightness galaxies which tend to be more 
compact (see also Sect. 4.5). 
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Fig. 2. The extracted spectra of all the galaxies for which we could obtain the velocity dispersion. 
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Fig. 3. Images of our galaxies from the CDFS field (3 panels on the left: HST+ACS+F435W, HST+ACS+F606W, 
HST+ACS+F850LP, from left to right) and from the 0055-269 field (rightmost panel, VLT+FORS 1 +Gunn z). Each box is 
3x3 arcsec. The numbers in the boxes are the last 3 digits of the object's name (see TableQ. 



Small sizes might also be due to the presence of an AGN 
in some of our galaxies, which could as well explain the large 
fraction of our objects showing [Oil] emission, particularly in 
the CDFS field. Therefore we have checked the presence of our 
sample objects in the catalogue of X-ray sources presented by 
Alexander et al. ( 2003 1, obtained from the deep Chandra obser- 
vations of the CDFS field by Giacconi et al. ( 120021 1. We found 
only one detection: galaxy CDFS_00467, which has a 0.5-8 
keV X-ray flux of 2.2 X 10~ 16 erg cm" 2 s" 1 . For an unobscured 
type 1 AGN with an X-ray/optical ratio of 1, the correspond- 
ing observed magnitude would be R~26, therefore contributing 
about 6% of the galaxy luminosity. However, since the source 
has a hardness ratio of «0.1, it is probably heavily obscured 
and much fainter in the /?-band. Therefore its influence on the 
profile should be completely negligible. 

4.2. The Fundamental Plane 

The FP is a powerful tool to study the evolution of early-type 
galaxies. Because of its dependence on galaxy luminosity, it is 
sensitive to recent star formation episodes. Therefore the FP 
could be useful to study the influence of the environment on 
the galaxy evolution as, for example, star formation induced by 
merging. 

We show in Figure [5] the position of our z~ 1 galaxies on 
the FP obtained by J0rgensen et al. ( 1996 hereafter JFK96) for 
the Coma galaxies at z - 0.0248 in the B-band. The FP at 



1 .5 
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1 .5 



Fig. 5. The fundamental plane as obtained by J0rgensen et al. 
( 1995a 1995b I for the Coma galaxies (triangles). Symbols are 
as in Fig.|4] The dashed line represents the best-fit plane to the 
Coma galaxies as seen edge-on. The small dots represent the 
GallCS model galaxies (see Section 5.2). 
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Table 2. The best fit parameters for the FP in the B-band 



Sample 


Ref. 








r 


Coma CI. 





1.20±0.08 


-0.83= 


=0.03 


-0.17±0.20 


0.88 <z< 1.3 


1 


0.88±0.16 


-0.63= 


t0.04 


0.46±0.38 



References. 0: JFK96, 1: this work. 

compared to the local one, keeps a remarkably small scatter, 
particularly for the more accurate CDFS data, and shows a clear 
offset, as already noticed previously (e.g. Franx et al. 2000 1, 
and likely also a different tilt, suggesting that the evolution of 
early-type galaxies depends on their size and/or mass and/or 
stellar population. 

In order to parametrize the offset and the rotation we have 
fitted a plane to the distribution of our data set, using the ex- 
pression (e.g. JFK96): 



logR e = alogcr +ySlog</) e + y, 



(1) 



where (I) e = L/2nR 2 , given the definition of the effective ra- 
dius, is the average surface luminosity within R e . We remind 
the reader that log(/) e = -0A((ji) e - C), where (I) e is in solar 
luminosities per square parsec, (ji) e is in magnitudes per square 
arcsecond and the constant C depends on the band, but not on 
the cosmology (C B = 26.982, C K = 24.982). 

The rms scatter of the logR e residuals from the best-fit FP 
is 0.11, while it becomes 0.17 when the Coma cluster param- 
eters are assumed in the fit, as done in Fig. [5] The best-fitting 
values of a, /3 and y, have been derived by minimizing the sum 
of the square distances from the plane, weighted by the obser- 
vational errors. They are listed in Table 13 together with those 
obtained for the Coma galaxies of JFK96 by means of the same 
procedure. The 90% confidence regions of the joint distribution 
of a and /3, for our z ~ 1 galaxies, and for the Coma cluster 
sample of JFK96 (Figure |6j are well separated, showing that 
no pair of (a, f3) exists that simultaneously provides an accept- 
able fit to the two data-sets. As expected, such a result is almost 
completely due to ACS data, which have significantly smaller 
errors for both R e and (I) e . However, we have verified that it 
still holds, if we use the values of R e and {I) e obtained from 
the de Vaucouleurs (n = 4) fit to the surface brightness profile. 
Thus our data show that at the 90% confidence level the FP ro- 
tates with redshift, even if larger samples are needed to better 
constrain its evolution. 

In order to study the evolution of the Ml Lb ratio in a way 
consistent with previous work, we assume that the dynamical 
mass of the galaxy is (Michard ll980> : 



M = 5 



cr 2 R, 



(2) 



This assumption is equivalent to assuming R 1/4 homol- 
ogy among early-type galaxies. We relax this assumption in 
Section 4.4. 

The evolution of the Ml Lb ratio obtained in this way is 
shown in Fig. [7] and Fig. [8] The spheroidal galaxies at z~l 
have a brighter fi-band luminosity for the same mass than low 
redshift galaxies by up to a factor of 8, particularly for the 
smaller mass galaxies. We note however that our sample goes 



-0.5 



-0.6 - 



-0.7 - 



-0.8 



-0.9 




Fig. 6. 1-sigma (solid) and 90% (dashed) confidence regions 
of the joint distribution of the FP parameters a and for our 
data-set of z~l galaxies, and for the Coma cluster sample of 
JFK96. 
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Fig. 7. Mass to light ratio vs. mass for our galaxies at z~l and 
for galaxies in Coma. Symbols are as in Fig.0] 



to smaller masses than the Coma sample of J0rgensen et al. 
l H995all99 5bi. The evolution in the M/L B ratio is similar to 
that of the massive galaxies in clusters for ~ 25% of our galax- 
ies, while it is stronger for the rest of our sample (see Fig. [8}. 

Actually, as Treu et al. (J2p05ai and van der Wei et al. 
(2005), we find that the evolution is different for galaxies 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Redshift 

Fig. 8. Evolution of the mass to light ratio for our galaxies. 
Symbols are as in Fig. |4] and the continuous line is the linear 
fit obtained by van Dokkum & Stanford (2003 1 for the evolu- 
tion of the mass to light ratio of cluster massive spheroids up to 
Z = 1.21. The other lines are obtained with simple stellar popu- 
lation models using the Chabrier ( 2003 1 IMF with a formation 
redshift of 5.0, 3.0, 1.5 and 1.0 (top to bottom). 



with different mass. If we fit a linear slope r\ to the evolu- 
tion of the Ai/L B ratio (Alog(M/L B ) = rjz), then we obtain 
77 = -0.59 + 0.03 for our 8 galaxies with M > 5 x 10 w Mq and 
77 = -0.86 ± 0.05 for our 9 galaxies with M < 5 x 1O 1O M . 
These two slopes are significantly different, showing that the 
evolution depends on the galaxy mass (see also Section 5.1). 
The slope for the more massive galaxies is close to the value 
77 = -0.46 ± 0.04, obtained by van Dokkum & Stanford 
(2003 1 for cluster massive galaxies (M > 10 n Alo). In fact 
we obtain 77 = -0.52 ± 0.04 for our 4 most massive galaxies 
(M>10 u Mq). 

Since our sample is selected in the K s band (approximately 
rest-frame /-band), the evolution of the Ml Lb ratio can hardly 
be due to selection effects, as could be in the case of selection 
in the observed optical bands, which for z~ 1 are very sensitive 
to star formation. 

4.3. Other Scaling Relations 

Figure[9]shows the relation between the velocity dispersion and 
the absolute magnitude in the B band - the so-called Faber- 
Jackson relation (Faber & Jackson ll976> . Here the trends are 
less clear than with the FP, also because even the low redshift 
objects show a large scatter. However z~l galaxies appear to 
be brighter than their z — counterparts with the same velocity 
dispersion. 




A 

1 .8 - 



1 .6 - 



-18 -19 -20 -21 -22 -23 

M B 

Fig. 9. The Faber-Jackson relation in the rest-frame B-band for 
our galaxies at z~l and for galaxies in Coma. Symbols are as 
in Fig. 01 small dots as in Fig. EI the dashed line is the relation 
obtained at z = by Forbes & Ponman ( 1999I I. 



A clearer trend is shown by the purely morphological rela- 
tion between the effective radius and the surface brightness, i.e. 
the Kormendy ( 1977 ) relation (Fig. llOl i. For a comparable size, 
Z~\ early-type galaxies are about 1-2 magnitudes brighter in 
surface brightness than the low redshift ones. This could in part 
be due to a selection effect: in fact Figure^|shows the limit for 
Mb < -20, which roughly corresponds to our apparent magni- 
tude limit of R ^ 24.4 for measuring the velocity dispersion. 
We also remark that the distribution of R e at z~l, although it 
partially overlaps the one for Coma, it is more concentrated to 
smaller sizes. 

Using the Zf-band absolute magnitudes computed by 
Fontana et al. (2004) by modest extrapolations of Bruzual & 
Chariot ( 2003 ) models fitted to the spectral energy distribution 
of our galaxies, we have also obtained the Faber-Jackson rela- 
tion in the A'-band (Fig. II li . Most of our galaxies are brighter 
also in the /T-band than their local counterparts with the same 
velocity dispersion, indicating some evolution of the M/Lg ra- 
tio. 

4.4. Masses 

We return to the determination of dynamical masses by relax- 
ing the assumption of homology, which was used in equation 
[2] In fact we know that our galaxies are not homologous, since 
they have different Sersic indices n, but seem to comply with 
a "weak homology", as n correlates with the luminosity (see 
Fig-El*- Bertin, Ciotti & Del Principe COO! BCD02 hereafter) 
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Fig. 10. The Kormendy relation in the rest-frame B-band for 
our galaxies at z~ 1 and for galaxies in Coma. Symbols are as in 
Fig.0]and small dots as in Fig. [5] The dashed line corresponds 
to M B = -20. 

have expressed the dynamical mass of a spheroidal galaxy in 
the form: 



Mk v = K v - 



cr 2 R„ 



(3) 



b 2.0 




-23 



-24 



-25 



-26 



-27 



-28 



Fig. 11. The Faber-Jackson relation in the rest-frame A'-band 
for our galaxies at z~l. Symbols are as in Fig.[4l the dashed 
line is the local relation obtained for Coma cluster galaxies by 
Pahre, Djorgovski and de Carvalho ( 1998 ). 



Fig. 12. The comparison of dynamical and stellar masses for 
our sample at z~l- Symbols are as in Fig.|4] 



where K v is a virial coefficient which takes into account the 
specific density distribution of both luminous and dark matter 
(DM), the specific star orbit distribution, and projection effects 
(see also Lanzoni & Ciotti 2003 1. By assuming spherical sym- 
metry, global isotropy of the velocity dispersion tensor, and ab- 
sence of DM, or a DM distribution which exactly parallels the 
stellar one, for a Sersic profile, Ky only depends on the Sersic 
index n and on the aperture within which the velocity disper- 
sion cr has been measured. For instance, the commonly used 
factor 5 of equation|2]is appropriate for a de Vaucouleurs R 1 ^ 
profile and for aperture radii of the order of R e /10. An accurate 
fit of Ky for 1 < n < 10 and small apertures (R e /&) can be 
found in BCD02. We have recomputed it for the 1.7 kpc aper- 
ture to which we correct our velocity dispersion measurements. 
The resulting values of Ky and of Mk v are listed in Tabled 
We remark that this dynamical mass is larger by a factor of up 
to 1 .55 than that estimated using equation|5] We emphasize the 
existence of an elliptical galaxy, i.e. CDFS_00633, at z=1.0963 
with red colors, a R 1 ^ 4 profile, no line emission, probably very 
low extinction, and a dynamical mass of 3.5 x 10 n /V(©. 

We have compared our dynamical masses obtained with 
equation [5] with the stellar masses obtained by Fontana et 
al. (2004 1 for the K20 galaxies using their Best Fit model, 
where the spectrum best fitting the complete U BVRIzJK s mul- 
ticolor photometry for Salpeter IMF is used. However since the 
Salpeter IMF is known to be inadequate (see e.g. Renzini 2005 
and Bruzual & Chariot 2003 1 and all empirical determinations 
of the IMF indicate that its slope flattens below ~ 0.5 Mq 
(Kroupa l200T1 Gould et al. ll996l Zoccali et al. l2000l . we have 
divided the stellar masses of Fontana et al. (2004i by a factor 
of 1.72, which is the age-averaged correction for a Chabrier 
(2003 1 IMF. The comparison is shown in Fig. 1121 
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The agreement between dynamical and stellar masses is 
good for the high mass galaxies, but most of our lower mass 
galaxies have stellar masses larger than the dynamical ones, 
with the maximum difference reaching a factor of four. Similar 
results have been recently found for lower redshift galaxies by 
Drory, Bender & Hopp (2004), who however attribute the dis- 
crepancy to the limited reliability of the SDSS velocity disper- 
sion measurements for low masses galaxies. We also remark 
that, given our large spectroscopic aperture (in units of R e ), 
in case of substantial rotational support, the adopted normal- 
ization of the observed velocity dispersion (Sect. 3.2) might 
be incorrect, and the dynamical masses of the stellar compo- 
nent underestimated (Riciputi et al. |2005 1. This might be the 
case particularly for the lower mass objects, that are also those 
showing a faster evolution in the Ai/L B ratio (see Section 5.1 
and Fig.1131. 

4.5. Comparison with previous results 

Van der Wei et al. ( 2005 1, who have also analysed the FP for 
galaxies in the CDFS, have two objects in common with us: 
their objects CDFS-18 and CDFS-19 are our CDFS.00633 
and CDFS_00571 respectively. Our results are in reasonable 
agreement with theirs. More specifically, considering ours vs. 
their results for CDFS.00633 and for CDFS_00571, the R e is 
6.65 + 0.12 vs. 4.17 ± 0.3 kpc and 6.81 ± 0.48 vs. 2.69 + 0.2 
kpc, the velocity dispersion is 260 + 23 vs. 324 ± 32 km s _1 
and 182 ± 21 vs. 229 + 35 km s" 1 , the K magnitude is 16.94 
vs. 16.98 and 17.54 vs. 17.53, the U - B color is 0.54 vs. 0.32 
and 0.26 vs. -0.32, the derived dynamical mass in I0 n Mq is 
3.5 ±0.7 vs. 5.1 ± 1.0 and 1.7 + 0.4 vs. 1.7 + 0.5 (we have used 
here the masses derived with equation|2] which is also assumed 
by van der Wei et al. l2005> . 

This comparison suggests that our effective radii are some- 
what larger than theirs, particularly for CDFS_00571. Our ve- 
locity dispersions and dynamical masses are comparable with 
theirs. Surprisingly they have much bluer colors than us, partic- 
ularly for CDFS_00571, for which their value seems too blue, 
and in fact it is not plotted in their Fig. 8. 

With respect to the morphological parameters obtained for 
the K20 galaxies in the CDFS field by Cassata et al. (I2005i . we 
have refined those previous values on a much smaller sample of 
objects by interactive optimization of fits of individual galaxies, 
masking out nearby objects and adjusting the local sky level. In 
general Cassata et al. (|2005 1 also find a decrease of the galaxy 
size with redshift, particularly for the elliptical galaxies. 

Also Daddi et al. ( 2005 1 have obtained morphological pa- 
rameters for 7 luminous early-type galaxies with 1 .4 < z < 2.5 
using Sersic profile fits with GALFIT on very deep HST+ACS 
images. They obtain small effective radii, even slightly smaller 
than ours, and a wide range of Sersic indices: 1 < n < 10. 

Finally, in comparing our relatively small sizes to the sizes 
obtained by others it is important to note that our spectroscopic 
galaxy type selection favours compact objects, while a morpho- 
logical type selection, such as e.g. that of Treu et al. ( 2005a|and 
2005b I, favours brighter and therefore probably larger objects. 



5. Discussion 

5.1. Spheroids Evolution 

The observed offset of the z~l FP from the local one is the 
result of the evolution of one or more physical properties of 
spheroids. This could be either the M/L ratio, as in the clas- 
sical analysis, or the size R e , or a combination of these two. 
Some evolution in the luminosity is actually expected as con- 
stituent stars grow older, even in the most passive scenario (see 
the models in Fig. |8j. However there is also evidence of size 
evolution both from our data and from others (see Section 4.1). 
It is however not clear what is the exact relative weight of these 
two contributions to the evolution of spheroids, and how does 
an individual object evolve in size. 

The observed rotation of the FP implies that the evolution 
of spheroids is not the same for all of them, whatever is the 
evolving parameter. Some have argued that the evolution could 
be mass dependent in the sense that lower mass galaxies evolve 
later than the massive ones (Fontana et al. 2004 Pozzetti et al. 
2003 1. Very recently Treu et al. (2005a) and van der Wei et al. 
(2005) have analysed the results of kinematic studies of field 
early-type galaxies with 0.2 < z < 1 .2 in term of a mass de- 
pendent evolution. In Section 4.2 we have already shown that 
the Ml Lb evolution of our early-type galaxies depends on their 
dynamical mass. This is clearer in Figure ^] where the differ- 
ence between the M/Lb evolution of our galaxies and the fit 
derived by van Dokkum & Stanford ( 2003 ) for the cluster mas- 
sive galaxies (6Alog(M/L B ) = Alog(M/L B ) +0.46z) is plotted 
as a function of mass. The dependence of the differential M/Lb 
evolution on mass appears to have a threshold around 10 u Mq, 
in the sense that more massive galaxies evolve in a way similar 
to that of cluster massive galaxies, while less massive objects 
evolve progressively faster and their M/Lb ratio is lower than 
that of cluster massive galaxies by a factor of up to about 3 for 
galaxy masses of the order of 10 1() A1o. We emphasize that the 
differential Ml Lb evolution, clearly visible in Figure ^]and in 
similar figures by Treu et al. (2005a) and by van der Wei et al. 
(2005), implies that the FP must rotate with redshift (Renzini 
& Ciotti fT^93l . 

The differential M/Lb evolution which we observe is sim- 
ilar to the so-called "downsizing" (Cowie et al. 1996 see also 
Kodama et al. 120041 Treu et al. 2005aJ, i.e. the later and/or 
longer lasting formation of lower mass galaxies. The reason for 
such later evolution of less massive spheroids could be found 
in a change with mass of either their stellar populations (IMF, 
age, metallicity), or of some structural/dynamical parameter, 
such as the DM distribution, the density profile, the degree of 
anisotropy and the partial rotational support, or a combination 
of these. However the only structural/dynamical changes possi- 
ble during the Universe lifetime are those induced by merging. 

A testable hypothesis is that the lower mass-to-light ra- 
tio of the less massive galaxies is due to recent star formation 
or younger ages increasing their B-band luminosity. We have 
therefore looked at the dependence of dAlogiM/Ls) on star 
formation indicators as the [Oil] equivalent width and the rest- 
frame U - B color. The dependence on the [Oil] equivalent 
width (Fig. [21 does not show any correlation, indicating that 



S. di Serego Alighieri et al.: The evolution of early-type galaxies at z~l. 



11 



0.4 



0.2 - 



0.0 



-0.2 



-0.4 - 



-0.6 



-0.8 




■.4 



0.2 



0.0 



-0.2 



-0.4 



-0.6 



-0.8 



9.5 



10.0 



10.5 11.0 
og Jl [M-a] 



1 1 .5 



4 6 
EW[0II][A] 



Fig. 13. The differential evolution of the Ml Lb ratio for each 
individual galaxy of our sample with respect to that of old clus- 
ter massive galaxies at the same redshift (see text). Symbols are 
as in Fig.|4] 



Fig. 14. The dependence of the differential evolution in Ml Lb 
on the rest-frame equivalent width of the [Oil] 3727 emission 
line for our spheroidal galaxies at z~ 1 . Symbols are as in Fig.0] 



the lower Ml Lb ratio is not due to on-going star formation ac- 
tivity. This is actually consistent with the fact that our galaxies 
have been selected as having early-type spectra and that the de- 
tected [Oil] line fluxes correspond to a rather low star forma- 
tion rate, ranging between 0.10 and 0.18 Mq yr _1 , using the 
conversion of Kewley, Geller & Jansen (2004) and assuming 
no extinction. On the other hand a correlation exists between 
the differential evolution in Ml Lb ratio and the U - B color 
(Fig.ll5>. with galaxies with lower Ml Lb ratio being bluer, in- 
dicating that indeed star formation activity has been proceeding 
until a rather recent past in the lower mass galaxies. 

As anticipated in Sect. 4.4, the combined effect of substan- 
tial rotational support and large spectroscopic apertures in the 
low-mass galaxies (e.g., Davies et al. 1983), could produce a 
trend similar to the differential evolution. In fact, aperture ve- 
locity dispersion of isotropic rotators, measured at or beyond 
R e , can be significantly lower than the same quantity mea- 
sured in the central region. Thus, the empirical correction of 
J0rgensen et al. i 1995b! see Sect. 3.2), that has been locally 
determined from velocity dispersion profiles obtained through 
small spectroscopic apertures (where the effect of rotation is 
negligible), might not be appropriate in our case (large aper- 
tures) if rotation is important. Unfortunately, a quantitative es- 
timate of such effect is not easy; however, a detailed analysis of 
rotational and flattening effects on aperture velocity dispersion 
of axysimmetric galaxy models shows that the variation is of 
the order of ss 30% for apertures of ss R e (Riciputi et al. 120051 
it increases up to x 80% when the whole galaxy image is taken 
into account, as assumed by Bender, Burstein & Faber 1992), 
much less than our estimated differential evolution. Note how- 
ever that the correction could be larger if a rotating disk is also 
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Fig. 15. The dependence of the differential evolution in Ml Lb 
on the rest-frame U-B color. Symbols are as in Fig. [4] 



present (see CDFS_00571 in Fig . EJ . In order to check the im- 
portance of this effect on our data, we have looked for signs 
of rotation on the two-dimensional spectra of our galaxies with 
line emission. Unfortunately the emission lines are too faint 
and compact to possibly show any sign of rotation at the S/N 
ratio of our spectra. 
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Finally, we remark that data in Fig. [O] were obtained un- 
der the assumption of strong R 1 ^ 4 homology (Eq. [2j, and so 
it is natural to ask how the same figure is modified by adopt- 
ing instead Eq.[3] In fact we know that our fainter galaxies have 
smaller Sersic indices (see Fig. 0}, and, since low n correspond 
to high values of Ky, this may reduce the predicted differential 
evolution. However the maximum effect due to non-homology 
is of the order of 55% (see Sect. 4.4), and therefore, as in case 
of rotation, much smaller than the required factor of 3. 

Of course, the two effects above could combine and reduce 
the amount of the required intrinsic differential evolution, but it 
seems to be impossible to blame just systematic effects for the 
behavior of data in Fig. ^] for which younger mean ages for 
low mass galaxies probably have a major effect (Fig.[T5l. 

5.2. Comparison with Hierarchical Models Predictions 

One of the currently open questions of observational cosmol- 
ogy is whether the observed evolution of early-type galaxies 
is consistent with the predictions of the hierarchical theories of 
galaxy formation. Our sample of early-type galaxies, being se- 
lected in the K band, which is particularly sensitive to galaxy 
mass, and being at the highest redshifts for which kinematic 
information has been obtained so far, is particularly well suited 
for this check. 

If the early-type galaxies form by merging, as predicted by 
the hierarchical models, then we should see a different evolu- 
tion in the cluster and in the field environment, where the merg- 
ing probabilities can be different. We fail to observe such effect, 
since our most massive field spheroids evolve at a rate indis- 
tinguishable from that of massive spheroids in clusters, and in 
the lower mass range there appear to exist early-type galaxies 
in clusters which evolve at a faster rate, comparable to that of 
our less massive field galaxies. Examples of this are our two 
galaxies at z = 0.67, which definitely evolve faster than cluster 
massive galaxies, although they are also likely in clusters. We 
note that also the two galaxies with M < K^'Al© in fig. 6a of 
van Dokkum & Stanford (2003 1 and the two lowest mass E+A 
galaxies observed by van Dokkum et al. i 19981 1 in a cluster at 
Z = 0.83 appear to have a lower Ai/L ratio than the massive 
galaxies in the same cluster. We do not exclude that the evo- 
lution of spheroids may also depend on the environment, but 
to test this hypothesis it is necessary to study the FP on well 
selected and complete samples of cluster early-type galaxies, 
overcoming the obvious observational effect, that usually pre- 
vents less luminous galaxies from finding empty slits in multi- 
object spectroscopic observations of high density cluster fields. 

If, as predicted by some versions of the hierarchical merg- 
ing model, the stellar populations of more massive galaxies 
are appreciably younger than those of smaller galaxies, the FP 
should become flatter with increasing redshift (e.g., Renzini 
1999 1. This is contrary to the observed steepening of the FP, 
and suggests that, if spheroids form by merging, the building 
blocks of the most massive objects at z ~ 1 should have formed 
and assembled well before those constituting the lower mass 
galaxies at the same redshift. 



To further check the ability of the hierarchical scenario 
in describing how galaxies form and evolve in the Universe, 
we have compared our results with the predictions of GallCS 
(Hatton et al. 2003 1, a hybrid (N-body and semi-analytic) hi- 
erarchical model of galaxy formation, that matches rather well 
the redshift distribution of the K20 survey (Blaizot et al. 2005 1. 
To perform a proper comparison with the observations, we 
have selected all field elliptical and SO galaxies with K s < 20 
and R < 24.4 in 3 timesteps, corresponding to redshifts z = 
0.87, 1.01 and 1.23. Although semi-analitic models are unable, 
by construction, to resolve the internal structure and dynamics 
of galaxies, spheroidal components in GallCs are approximated 
as Hernquist ( 1990 1 spheres, characterized by a half-mass ra- 
dius r/, and a velocity dispersion at r/, (see Hatton et al. {2003 1 
for details). We have therefore used the Hernquist model to 
correct these values and derive the corresponding effective ra- 
dius R & and the one-dimensional velocity dispersion within a 
1.19 h" 1 kpc aperture, consistently with what done for the ob- 
servations. The resulting values and the corresponding effective 
surface brightness have been used to build the simulated scal- 
ing relations, that are shown in Figures [5] |9] and ^3 The pre- 
dicted FP appears to be similar to the observed one, even if its 
scatter is larger than observed. A good agreement is found for 
the Kormendy relation, while the discrepancies at the faint end 
of the Faber-Jackson relation suggest excessively large veloc- 
ity dispersions and/or underestimated /^-luminosities for sim- 
ulated objects with Mb $ -21. The model stellar masses, in- 
stead, are found to be in the correct range, between 10 10 and 
2 x 1O U A10. If one considers that the definition of morpho- 
logical types in the model (based on the bulge-to-disc lumi- 
nosity ratio in the B band) is different from the one adopted 
in the observations, and that structural and dynamical prop- 
erties of galaxies can be only roughly described in this kind 
of modelling, the present comparison can be considered quite 
satisfactory. However, a major problem exists concerning the 
colours of simulated galaxies, which are too blue compared 
to the observations: for the selected GallCS sample, the rest- 
frame B - K ranges between 1.9 and 4.6, with median at ~ 3, 
while all observed galaxies have B — K > 3.2. In addition, ex- 
cessively large dust extinctions are needed to obtain the reddest 
colours in the model galaxies. This is a well known, but still un- 
solved, problem common to all semi-analytic and hybrid mod- 
els of this kind (see, e.g., Firth et al. 2002, Benson et al. 2002 
Pozzetti et al. 2003, Somerville et al. 2004 1, suggesting that a 
re-examination of how star formation proceeds in galaxies is 
required, particularly at high redshift. 

6. Conclusions 

We have obtained photometric (total magnitudes in various 
bands), morphological (effective radius, average surface bright- 
ness within the effective radius, Sersic index) and kinematic 
(velocity dispersion) parameters for an almost complete sample 
of 15 field early-type galaxies selected from the K20 survey in 
the redshift range 0.88 < z < 1.30. Apart from completeness, 
our sample has the advantage of being selected in the ^T-band, 
of having the galaxy type assigned spectroscopically, and of 
covering two separate fields. 
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From these data we obtain the following results: 

1. The distribution of effective radii is shifted towards 
smaller sizes compared to the local one. 

2. Sersic indices span a range between 0.5 and 4, and cor- 
relate with B-band luminosity. 

3. The resulting FP in the rest-frame B-band at z~l is offset 
from, and steeper than, the local one, but keeps a remarkably 
small scatter. 

4. Under the assumption of homology, and at fixed R e , the 
evolution of the dynamical Ml Lb ratio shows a dependence 
on galaxy mass: for our more massive objects, it is similar to 
that of cluster massive galaxies, while it is faster for the less 
massive ones, with a threshold around 10" Mq. 

5. For a given velocity dispersion or size some of the z~l 
galaxies are considerably brighter than local ones, both in the 
rest-frame B- and /f-band. 

6. Dynamical masses (properly obtained by taking into ac- 
count the observed Sersic index) are consistent with stellar 
masses (derived from model fits to the spectral energy distri- 
bution) in the high mass range, while they are smaller in the 
low mass range. 

7. In about half of our z~l early-type galaxies we have de- 
tected [Oil] line emission, corresponding to a low star forma- 
tion rate ranging from 0.10 to 0.18 Mq yr~ . 

We give the following interpretation: 

1. Spheroidals evolve in M/L ratio, and possibly also in 
size. 

2. The rate of evolution over the last 2/3 of the Universe 
lifetime depends on the galaxy mass, low mass galaxies evolv- 
ing faster. 

3. A fraction of the detected differential evolution might 
be due to a stronger rotational support in low mass galaxies 
with respect to the high mass ones, and to the fact that obser- 
vational data are routinely interpreted under the assumption of 
strong homology. However the correlation with galaxy color 
indicates that a significant fraction of the detected differential 
evolution is indeed intrinsic, i.e. due to more recent star forma- 
tion episodes (the so called "downsizing"). 

4. We fail to detect differences in the evolution with the 
environment, and the rotation of the FP goes in the opposite 
sense than predicted by the hierarchical merging models. 

5. The GallCS hierarchical galaxy formation model is able 
to produce galaxies which are in the correct mass and size range 
and reproduce the various scaling relations, although with a 
larger scatter for the FP and with unrealistic colours and dust 
extinctions. 

Since neither homology (constant Ky) nor the other ex- 
treme case of weak homology (constant mass-to-light ratio, 
BCD02) appear to be holding for our z~ 1 spheroids, we are 
astonished that the FP keeps such a small scatter, while shift- 
ing and rotating as the Universe evolves, thus posing additional 
and more stringent challenge to our understanding of galaxy 
formation and evolution. 
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